A blackbody is not a blackbox 
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We discuss carefully the blackbody approximation, stressing what it is (a limit case of radiative 
transfer), and what it is not (the assumption that the body is perfectly absorbing, i.e. black). 
Furthermore, we derive the Planck spectrum without enclosing the field in a box, as is done in most 
textbooks. Athough convenient, this trick conceals the nature of the idealization expressed in the 
concept of a blackbody: first, the most obvious examples of approximate blackbodies, stars, are 
definitely not enclosed in boxes; second, the Planck spectrum is continuous, while the stationary 
modes of radiation in a box are discrete. Our derivation, although technically less elementary, is 
conceptually more consistent, and brings the opportunity to introduce to students the important 
concept of local density of states, via the resolvent formalism. 



"No colors anymore, I want them to turn black 
[...] I wanna see the sun blotted out from the sky 
I wanna see it painted, painted, painted, painted black 

Yeah!" 

The Rolling Stones 



I. INTRODUCTION 

Planck's law for blackbody radiation is celebrated as a 
landmark in the history of physics, being the first physi- 
cal law conflicting patently with classical mechanics. De- 
rived in 1900, it is usually regarded as the igniter of the 
century of quantum mechanics. But its relevance goes 
beyond quantum mechanics. Blackbody radiation is a 
striking example of universality in statistical mechanics 
- the light radiated by a blackbody does not depend on 
its constitution, only on its temperature -, a topic which 
came to the foreground with the later studies of criti- 
cal phenomena. More recently, Planck's distribution has 
appeared to fit with an astounding accuracy the cosmic 
microwave background (deviations are at most 50 per 
millioni), thus opening the era of precision cosmology. 

Because of its seminal importance, one feels that 
the derivation of Planck's law presented to students 
should be as lucid as possible. In most textbook 
presentations^r— , however, one apparently innocuous step 
undermines the understanding of its applicability: the 
radiation field is assumed to be enclosed in a box. This 
prompts students to think that the Planck spectrum is 
the one radiated by a black box. This is true, of course, 
but also very misleading. If the blackbody radiation is 
really the blackbox radiation, why should stars, which 
are not enclosed in boxes, have a Planckian spectrum? 
Why should hot metal bars, bulb filaments, incandescent 
lava or the background of the universe have a Planckian 
spectrum? 

Since the physics of thermal radiation generally does 
not involve a box, we feel that the discussion of Planck's 
law should not either. In this note, we propose such a 
discussion. We introduce the blackbody approximation 



in the framework of radiative transfer theory, and derive 
the Planck spectrum using the notion of local density of 
states. We do not claim originality in the method used, 
the resolvent formalism, which is standard in condensed 
matter physics and scattering theory. Rather, our aim is 
to promote an approach we believe to be pedagogically 
more transparent. (See also Refi^ for a more advanced 
discussion of Planck's law from the many body theory 
perspective.) 

The paper is organized as follows. In the next section, 
we discuss the blackbody approximation, and its relation 
with Kirchhoff's law. In section III, we define the local 
density of states, and derive Planck's law without the 
box, using the resolvent formalism. Section IV presents 
our conclusion. 



II. WHAT IS A BLACKBODY? 

A. The common definition 

A blackbody is usually defined as "a body which com- 
pletely absorbs all radiation incident on it" . Landau and 
Lifschitz^ add 

Such a body can be realised in the form 
of a cavity with highly absorbing internal 
walls and a small aperture. Any ray enter- 
ing through the aperture can return to it and 
leave the cavity only after repeated reflection 
from the walls of the cavity. When the aper- 
ture is sufficiently small, therefore, the cavity 
will absorb practically all the radiation inci- 
dent on the aperture, and so the surface of 
the aperture will be a black body. 

Clearly, this picture makes perfect experimental sense. 
Theoretically, on the other hand, the nature of the ide- 
alization it is meant to express remains clouded. If a 
blackbody emits light, why is it called "black"? And if 
it is an ideal emitter, why is it defined as an ideal "ab- 
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sorber"? What is the role of the cavity? More impor- 
tantly, what would an imperfect blackbody be? A cavity 
which does not absorb all incident radiation? One with 
a bigger aperture? 

B. A continuous spectrum 

Many objects emit light like imperfect blackbodies. 
The most prominent among them is, of course, our Sun 
(see Fig. 1). Its spectrum, first measured in the early 
nineteenth century by Wollaston and Frauhofer, is well 
approximated by that of a blackbody of temperature 
between 5500 and 6000 K. It does displays deviations 
from it, the so-called Fraunhofer lines (corresponding to 
the absorption of certain frequencies by the solar atmo- 
sphere), but it is indisputable that, athough very dissim- 
ilar from an absorbing cavity with a tiny hole, the Sun 
qualifies as an imperfect blackbody. On the other hand, 
the spectroscopy of monoatomic gases shows clearly that 
their emission spectrum is not Planckian. They emit only 
certain discrete frequencies (Fig. 2), determined by the 
electronic structure of the atomsi^i 

This observation hints at what a blackbody really is: a 
body with a rich energy spectrum, capable of exciting all 
frequencies of light by thermalization.^^ 

From this perspective, the box definition appears para- 
doxical. As is well-known, a closed box selects certains 
light frequencies, through the condition uji = mrc/Li, 
where Li is the dimension of the box in the z-direction. 
Thus, instead of permitting a wide range of thermally 
excited frequencies, the box restricts the emission spec- 
trum, even making it discrete. Of course, one could argue 
that the volume of the box can be made arbitrarily large, 
and therefore that this quantization of frequencies is not 
physically relevant. But this is precisely our point: as far 
as the frequencies of light are concerned, the box is not 
physically relevant. 

C. Kirchhoff's law 

In its standard definition, a blackbody is one that "ab- 
sorbs all incident light" - a black body. This fact alone 
should disturb the mindful student: how can a body be 
black, and yet emit a colourful spectrum of thermal light? 
In any case, what should the absorptive power of a body 
have anything to do with the caracteristics of its thermal 
emission? 

The answer lies in an experimental observation which 
played a key role in the nineteenth century developments 
which led to Planck's successful analysis of thermal radi- 
ation, and which is too seldom mentioned in undergrad- 
uate discussions of thermal radiation - Kirchoff's lawJ^ 

The radiative properties of a body are characterized by 
its emissivity and absorptivity (and scattering, which can 
usually be neglected). These can be defined by the fol- 
lowing schematic model for the propagation of radiation 
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FIG. 1: The Sun is an imperfect blackbody, monoatomic gases 
are not (top: hydrogen; bottom: iron). 

within a mediumfii As the (monochromatic) beam trav- 
els trough the medium, the variation of its energy density 
u{l,Lj) receives two contributions: a positive one, corre- 
sponding to emission, and a negative one, corresponding 
to absorption: 

dvL 

— {l,uj)=£{uj)-a{uj)u{l,uj). (1) 

The coefficients £{uj) and a{uj) are the emissivity and ab- 
sorptivity of the body. Now, Kirchoff's law states that, 
although e(w) and a(w) largely depend on the constitu- 
tion of the material, at thermal equilibrium, their ratio 
Jt{uj) = e{uj) / a{uj) is universal; it depends on temper- 
ature and frequency only. A good absorber (a(w) large) 
at a certain frequency is also a good emitter {e{uj) large) 
at that frequency, and vice versa. 

At this point, our mindful student's worries should al- 
ready be eased: if Kirchhoff's law is right, then black- 
bodies, which are by definition excellent absorbers, must 
also be excellent emitters. But further reflection should 
reveal a caveat in this line of thought: emissivity and 
absorptivity usually depend on the actual material used, 
while the blackbody radiation it emits does not. Why is 
that? 

Some insight into this question is provided by the fol- 
lowing consideration, due to EinsteinJ^ The interaction 
between matter and radiation boils down to transitions 
between energy levels: given two levels a and b, with 
respective energies Ea < Ef,, emission corresponds to 
an upgoing transition a ^ b, while absorption corre- 
spond to a downgoing transition b a. The rate of 
these transitions, Fa^f, and Ff,_).(j, is what controls at the 
microscopic level the absorptivity and emissivity of the 
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body. Now, the condition of thermal equilibrium fixes 
the probabilites of each levels, through the Gibbs distri- 

bution pa^b oc e ''bt _ But, and this is the key point, it 
has no bearing on the transition rates Ta^b and Tb~^a 
themselves, but only on their ratio. Indeed, in this mi- 
croscopic perspective, thermal equilibrium translates into 
the condition of detailed balance, according to which the 
probability flux between microstates cancel exactly: 



This leads to 



Pa^a- 



■ Pb^, 



b^a- 



(2) 



(3) 



absorbers (a(a;) large), but also for objects involving long 
optical paths {I large). 

This is precisely what a closed cavity does: it provides 
the conditions for this convergence process to take place, 
whatever the intrinsic properties e and a of the material 
within the cavity. The box, in other words, is a useful ex- 
pedient for the actual production of thermal radiation in 
the laboratory when the hot body is not black: it permits 
to reach the optically thick limit even if the absorptivity 
is low, because it generates very long paths I withing the 
material. But the presence of a box is of course not a 
necessary condition for the thermalization of light: stars 
too are optically thick - and this is why their spectrum 
looks Planckian. 



in which the RHS is a function of the transition frequency 
LOab and temperature only. Kirchhoff's law is a conse- 
quence of this constraint on the transition rates imposed 
by detailed balancei^ 

Combining the phenomenological radiative transfer 
equation ([T]) with Einstein's microscopic model, we un- 
derstand that thermal equilibrium of the material source, 
through the condition of detailed balance, constrain the 
ratio of emission to absorption - Kirchhoff's law - but 
not their respective rates independently: this is why dif- 
ferent materials have different thermal emission and ab- 
sorption properties. It is only under a further assumption 
that a universal function describing blackbody emission 
- Planck's function - can be obtained. 

What is this further assumption? Is it that the body 
"absorb all light incident on it" , as in the standard defini- 
tion? In other words, that the body be perfectly absorb- 
ing ('black') on the whole spectrum? No! Such a material 
does not exist.— That is, the condition of 'blackness' en- 
tering the standard definition of a blackbody cannot be 
the idealization underlying blackbody radiation. Then 
what is it? 



D. Optical thickness 

The answer should be obvious by now: the additional 
condition is that the radiation field itself should be at 
thermal equilibrium. This means that the random pro- 
cesses of emission and absorption of light by the hot body 
should have reached their stationary state. Solving ([T]) 
with the initial condition u{l — — 0, we find 



= n _ 



cj) = (1 — e 



)JtH 



(4) 



In particular, we see that if a{uj)l ^ 1, the so-called 'op- 
tically thick' limit, the observed energy density is given 
by Jt{^jj), which is nothing but the Planck spectrum. In 
other words, the coefficient a{uj) measures the rate of 
convergence of the emitted light to its equilibrium value 
Jt{'^) ~ the blackbody approximation is just the condi- 
tion that a{uj)l ^ 1. This is achieved not only for good 



III. PLANCK'S LAW WITHOUT THE BOX 

A. Partition of energy 

From the previous discussion, we know that the equilib- 
rium spectral energy density coincides with the function 
Jt in Kirchoff 's law: the energy contained in thermal ra- 
diation within an infinitesimal volume d^x about a point 
X in the frequency range duj is 



dE{x,uj) = JT{x,Lj)d^xduj 



(5) 



How can one compute this function? The answers to 
this question is obvious when one recalls that the total 
energy of the field is distributed over its modes, and that 
modes with the same frequency lo are degenerate, i.e. 
store the same energy. Hence: 

dE{x,uj) — (energy of each w-mode) 

X (occupation number of each w-mode) 
X (number of w-modes accessible at x) (6) 

From quantum mechanics, we know that the energy 
stored in a mode with frequency uj is E{uj) = huj. More- 
over, the occupation number ^^(a;) of each such mode at 
thermal equilibrium is given by the Bose-Einstein distri- 
bution, 



1 



(7) 



It follows that the thermal energy density Jt{x,uj) is 
given by 



Jt{x, uj) — hu! 



1 



pix,uj), 



(8) 



where p{x, lo) denotes the density of modes with fre- 
quency UJ around the point x. Indeed, for an observer 
localized at x, certain modes might be inaccessible, or 
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only partly accessible: p{x, uj) is the space-resolved den- 
sity of modes, often called (using quantum mechanical 
parlance) the local density of states (LDOS). Of course, 
if the medium is homogeneous, this quantity does not 
depend on x. However, in more general situations (such 
those described at the end of the next section), it does. 

All in all, the Planck spectrum can be described as the 
outcome of the balance between the energy per mode 
and the LDOS, which tend to grow as uj gets large, 
and the Bose-Einstein distribution, which favors the low- 
frequency modes. This is illustrated in Fig. [21 




FIG. 2: Balance between the LDOS, the energy per mode and 
the Bose-Einstein distribution. 



B. A hint of spectral theory 

The reason why the evaluation of the LDOS is non- 
trivial mathematically is because the light frequencies 
span a continuous spectrum. Had it been discrete, the 
LDOS would just have been the degeneracy of w, under- 
stood as an eigenvalue of the wave operator. In the con- 
tinuous case, however, the notion of degeneracy is sub- 
tler, for it involves the definition of a measure on the 
spectrum. This is analogous to the ambiguity one faces 
when trying to extend a discrete sum to an integral: the 
appropriate weighting of the integration variable is not 
obvious anymore. The case of the standard measure on 
the sphere, which involves a factor sin0, is an obvious 
example of such subtleties. Of course, in this case, the 
appropriate measure is determined by the condition of ro- 
tational invariance. What determines the spectral mea- 
sure of an operator? 

A powerful tool to answer such a question is the 'resol- 
vent formalism']^ Given a (self-adjoint) operator A, the 
resolvent is the operator- valued function of the complex 
variable z defined by 



By construction, R{z) is analytic in C\a{A), the com- 
plement of the spectrum of A. The discrete part of (j{A) 
corresponds to isolated poles of R{z), while the contin- 
uous part of (y{A) generates a branch cut along the real 
axis. In other words, the spectrum of A is encoded in the 
analytic structure of R{z). 

Given any two vectors ip and (p, it follows from the 
spectral theorem that 



{^\R{z)ij) = 



E 



X-z 
.9(A) 



(10) 



X-z 



<Tc(A) 



p^,4i{X)dX 
X-z 



where dfj.^^^ is the spectral measure associated to ip and 
-0, and in the second equality the spectrum is decomposed 
into its discrete ad{A) and continuous ac{A) partsi^ The 
spectral density Pip,^,(A) is therefore defined as the den- 
sity of the continuous part of the spectral measure dfj,^^^ 
with respect to the Lebesgue measure dX. 

When ^ is a wave operator, and {ip, ip) = (x, y) 
are position (generalized) eigenvectors, the quantity 
G{x,y;z) = {y\R{z)x) is called in the physics literature 
the 'Green function'. The spectral decomposition of its 
diagonal elements G(x, x; z) reads 



G{x,x;z)^ ^ 



.9(A) 



A 



Px{X)dX 
X-z ■ 



(11) 



This expression provides the spectral density Px{^) with 
the following interpretation: the generalized eigenvalue 
A, when analyzed through a state localized at the point 
x, comes with a weight Px{X). This weight is the overlap 
of the density of A-eigenmodes with the local state |a;) - 
in other words, the LDOS p{x, A) of A. 

Moreover, the formula (|lip indicates a procedure to 
evaluate Px{X)- Indeed, just as the degeneracy g[X) of 
an isolated eigenvalue A can be read off as the residue 
of {y\R{z)x) at z = A, the LDOS is given by the dis- 
continuity of G{x, x; z) along the branch-cut singularity. 
This is the so-called Stieltjes-Perron inversion formula 
(Appendix A): 

p(a;,A) = -^limA,G(x,a;;A), (12) 

with 



R{z) = {A-z)- 



(9) 



A^G{x, x; A) = G{x, x; A -I- ie) — G{x, x; X — ie). (13) 



C. Resolvent of the Laplace operator 

Let us apply this resolvent formalism to the case of 
electromagnetism. For the sake of simplicity, we shall 
consider here the scalar Helmholtz equation ~ the two 
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independent polarizations of the field will be taken into 
account a posteriori, by multiplying the LDOS thus ob- 
tained by a factor of 2. 

The Helmholtz equation in vacuum for a monochro- 
matic wave 4>u, reads 

-A0„(x) = ^^^^(a:), (14) 

which is an eigenvalue equation for the Laplacian, with 
eigenvalue A = w^/c^. We introduce the resolvent 
(—A — z)~^, and evaluate the corresponding Green func- 
tion G{x, y; z) by Fourier transform 

G(-,2/;^)=y (27r)3 fc2 _ , ■ (15) 
D. The vacuum LDOS 




FIG. 3: LDOS versus frequency at different heights above a 
semi-infinite sample of aluminium, from Joulain et aliU The 
resonance around 10^^ rad/s corresponds to the excitation of 
surface-plasmon polaritons. 



The two integrals G{x, y; A ± ie) are readily evaluated 
by residue calculus (Appendix B), yielding 



lim A^G{x, x; A) 



iVX 



Using and the fact that dA = iLojc^dio, we obtain 
the electromagnetic vacuum LDOS (with a factor of 2 for 
the two polarizations): 



p(a;,w) 



(17) 



Note that, as expected, the translational invariance 
of the vacuum translates into the independence of the 
LDOS on the space point x. Plugging this value into ([8]) 
immediately yields Planck's law. 



(18) 



gKgJ 



E. Surface effects 

We conclude this section by an illustration of the 
strength of the LDOS method for a finer analysis of ther- 
mal radiation, in particular close to a metallic surface. 
At a distance to the hot body comparable to the thermal 
wavelength, the vacuum approximation breaks down, as 
the field reveals evanescent modes and polariton excita- 
tions. Such local excitations have been shown to modify 
significantly the LDOS, and hence the hot body spec- 
trum, notably by enhancing monochromaticity (Fig. 3), 
spatial coherence, and directivity^^"— This shows an- 
other facet of the blackbody approximation which is not 
usually emphasized: it is a far-field approximation. 



IV. CONCLUSION 

In this paper, we have tried to disentangle two aspects 
of the standard picture of a blackbody as a cavity with a 
small aperture. The first one relates to the efficiency of 
the thermalization of light through its almost everlasting 
interaction with the walls of the cavity: this is the one 
invoked in Landau and Lifschitz's definition. Although it 
does not make the conditions of thermalization of radia- 
tion explicit, the cavity picture is useful to demonstrate 
how a blackbody can be realized in the laboratory. The 
second aspect relates to the selection of certain frequen- 
cies, and the subsequent Fourier space mode-counting ar- 
gument. Unlike the former, this aspect does not merely 
serve illustrative purposes, but usually enters the actual 
derivation of the Planck spectrum. As such, it appears 
to students as an important feature of thermal radiation. 
We have argued that it is not, and that it actually con- 
tradicts a key feature of blackbody radiation, namely the 
fact that all wavelengths are emitted. 

Our approach, focused on local properties of the elec- 
tromagnetic field, allows to derive Planck's law without 
using a blackbox, and provides tools for subtler consider- 
ations, such as the near-field regime of thermal radiation. 
With this perspective, we believe that students are less 
likely to miss the point of the blackbody approximation, 
which expresses an idealization of the interactions be- 
tween the electromagnetic field and heated materials - 
and not of the materials themselves. 
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Appendix A: The Stieltjes-Perron formula 



Let a complex function / be defined as 

Pit) 

t - z 



fiz) = / dt 



(Al) 



where p is a continuous density on the open interval /. 
We consider the values = /(A ± ie), for A € /, and 
use the well-known identity (the Sokhatsky-Weierstrass 
theorem) 

lim , . = V{-^) ± 17:6(1 - A) (A2) 

where Vij^j) is the Cauchy principal value. It follows 
that 

lim (/+ - /") = / dt p{t)2iTrS{t - A) (A3) 
Ji 

= 2*^p(A), (A4) 
which is the Stieltjes-Perron inversion formula. 



Appendix B: Evaluation of G{x,y; A ± ie) 
Consider the two integrals 

f d^k Qik-{3:-y) 

G{x,y■,X±^e) = J (,^)3 fc2 _ ^ ■ (Bl) 

denoted for simplicity. The angular integration is 
straightforward, and yields 

The remaining integral is evaluated using the residue the- 
orem, considering as integration contours the standard 
half-circles 7±, enclosing the poles k± = VA ± ie respec- 
tively. We obtain 



g± ^ ^ i sm{k±\x-y\) 
47r |a; — y| ' 

The discontinuity across the cut is then given by 



(B3) 



lim(G+-G-) = f ""(/^'"7^'\ (B4) 
e^o^ ' 2t: \x-y\ ^ ' 

from which (|16p follows by setting y — x. 
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We disregard subtleties related to singular spectra. 



